Perturbative and non-perturbative parts of eigenstates and local spectral density of 
states: the Wigner band random matrix model 
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A generalization of Brillouin- Wigner perturbation theory is applied numerically to the Wigner 
Band Random Matrix model. The perturbation theory tells that a perturbed energy eigenstate 
can be divided into a perturbative part and a non-perturbative part with the perturbative part 
expressed as a perturbation expansion. Numerically it is found that such a division is important 
in understanding many properties of both eigenstates and the so-called local spectral density of 
states (LDOS). For the average shape of eigenstates, its central part is found to be composed of its 
non-perturbative part and a region of its perturbative part, which is close to the non-perturbative 
part. A relationship between the average shape of eigenstates and that of LDOS can be explained. 
Numerical results also show that the transition for the average shape of LDOS from the Breit-Wigner 
form to the semicircle form is related to a qualitative change in some properties of the perturbation 
expansion of the perturbative parts of eigenstates. The transition for the half-width of the LDOS 
from quadratic dependence to linear dependence on the perturbation strength is accompanied by 
a transition of a similar form for the average size of the non-perturbative parts of eigenstates. For 
both transitions, the same critical perturbation strength At has been found to play important roles. 
When perturbation strength is larger than A;,, the average shape of LDOS obeys an approximate 
scaling law. 

PACS number 05.45. +b 



I. INTRODUCTION 

The average shape of energy eigenfunctions (EFs), 
which characterizes the spreading of the eigenstates of a 
perturbed system over the eigenstates of an unperturbed 
system, is very important in a wide range of physical 
fields, from nuclear physics and atomic physics to con- 
densed matter physics. For example, in a recent approach 
for the description of isolated Fermi systems with finite 
number of particles, a type of "microcanonical" partition 
function has been introduced and expressed in terms of 
the average shape of eigenfunctions Q. However, up to 
now only some of the general features of the shape have 
been known clearly. For example, when perturbation is 
not strong, generally the shape has a Breit-Wigner form 
(Lorentzian distribution) |^-|^|, but in a larger range on 
the interaction strength, without any simple analytical 
expression known for the shape, only some phcnomcno- 
logical expressions have been suggested H [llj . 

Another quantity, the so-called local spectral density 
of states (LDOS), has also attracted lots of attention 
recently ]^,|^,|l2|-|ll| . This quantity, known in nuclear 
physics as the "strength function" , gives information 
about the "decay" of a specific unperturbed state into 
other states due to interaction. In particular, the width 
of the LDOS defines the effective "lifetime" of the unper- 
turbed state. For sufficiently weak coupling, the average 
shape of LDOS is usually close to the Breit-Wigner form 
with the half-width having a quadratic dependence on 



the interaction strength. On the other hand, when inter- 
action is strong, the shape is model dependent and the 
half-width is linear in the interaction strength. 

Numerically it is already known that for Hamiltonian 
matrices with band structure generally both the aver- 
age shape of EFs and that of LDOS can be divided 
into two parts: central parts and tails with exponen- 
tial (or faster) decay. However, an analytical definition 
for such a division has not been achieved yet. A possi- 
ble clue for this problem comes from a generalization of 
the so-called Brillouin- Wigner perturbation theory |l7f] 
introduced in Ref. [[l8| for studying long tails of EFs, 
which tells that analytically EFs can be divided into per- 
turbative and non-perturbative parts with perturbative 
parts being able to be expanded in a perturbation expan- 
sion by making use of the corresponding non-perturbative 
parts. The relationship between central parts and non- 
perturbative parts of EFs was not studied in Ref. Jl8| , 
while it is more important. 

The purpose of this paper is to study if the above 
mentioned generalization of Brillouin- Wigner perturba- 
tion theory (GBWPT) can provide deeper understanding 
for properties of EFs and LDOS, especially the relation- 
ship between the two divisions mentioned above for EFs. 
As a first step, here we study the so-called Wigner band 
random matrix (WBRM) model, which is one of the sim- 
plest models for the application of the GBWPT. 

The WBRM was first introduced by Wigner 40 years 
ago p3] for the description of complex quantum sys- 
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terns as nuclei. It is currently under close investigation 
]6|, ^2|Jl^j20| -p5t since it is believed to provide an adequate 
description also for some other complex systems (atoms, 
clusters, etc.) and as well as for dynamical conservative 
systems with few degrees of freedom, which are chaotic 
in the classical limit. Unlike the standard Band Ran- 
dom Matrices (see, for example, ]2^ , p7[ |), the theory of 
WBRM is not well developed. Numerically it is known 
that different averaging procedures for EFs give different 
results for their average shape. Specifically, when aver- 
aging is done with respect to centers of energy shell, for 
strong interaction central parts of averaged EFs are close 
to a form predicted for the LDOS by the semicircle law. 
On the other hand, for the case of averaging with respect 
to centroids of EFs, central parts of averaged EFs are 
obviously narrower than the corresponding LDOS when 
the so-called Wigner parameter is large and an ergodicity 
parameter is small 

Comparatively, more properties are known for the av- 
erage shape of the LDOS of the WBRM. For the tails, a 
corrected analytical expression has been derived in Ref. 
Q for weak perturbation and numerically found correct 
for even strong perturbation |q jll] , ^2[] . A more general 
analytical approach for the LDOS of the WBRM is given 
in Ref. {Hq] , where an expression for the LDOS is given in 
terms of a function satisfying an integral equation. In suf- 
ficiently weak and extremely strong perturbation cases, 
the expression gives the well known Breit- Wigner form 
and semicircle law, respectively. For the transition from 
the Breit- Wigner form to the semicircle form, although 
the expression also predicts correct results for the LDOS, 
the physical explanation for it is not so clear yet. There- 
fore, it is of interest to study if the GBWPT can throw 
new light on the above mentioned properties of the EFs 
and the LDOS of the WBRM. 

This paper has the following structure. For the sake 
of completeness, in section II we give a brief presentation 
of the GBWPT. Some predictions for properties of EFs 
of the WBRM are also given in the section. Numerical 
results for the size of non-perturbative parts of EFs are 
discussed in section III. It is shown that for the average 
shape of EFs the averaged non-perturbative part is in- 
deed related to the central part. A region predicted in 
section II, which is between the non-perturbative part 
and the (exponentially or faster decaying) tails of an EF, 
is also studied numerically. Section IV is devoted to a dis- 
cussion for the relation between the average shape of EFs 
and that of the LDOS. In section V, some properties of 
the average shape of LDOS, such as the transition from 
the Breit- Wigner form to the semicircle form and the 
dependence of its half-width on perturbation strength, 
are studied numerically and found to be closely related 
to properties of the average size of the non-perturbative 
parts of the corresponding EFs. Finally, conclusions are 
given in section VI. 



II. GENERALIZATION OF BRILLOUIN- WIGNER 
PERTURBATION THEORY 

For the sake of completeness, in this section we first 
give a brief presentation of the above mentioned gener- 
alization of Brillouin- Wigner perturbation theory (GB- 
WPT) pit]. Then, we give a brief discussion for some 
properties of the EFs of the WBRM. 

Generally, consider a Hamiltonian of the form 

H = H° + XV (1) 

where H° is an unperturbed Hamiltonian, V is a pertur- 
bation and the parameter A is for adjusting the strength 
of the perturbation. For the WBRM with dimension N 
and band width b, the Hamiltonian matrix considered 
here is chosen of the form 

H i:j = (H° + XV)ij = E ^ + Xvij (2) 

where 

E°i=i (i = l,---,N) (3) 

are eigenenergies of the eigenstates of H° labeled by \i). 
Off-diagonal matrix elements mj = Vji are random num- 
bers with Gaussian distribution for 1 < \i — j\ < b 
(( v ij) — an d (v?j) = 1) and are zero otherwise. Here 
basis states have been chosen as the eigenstates \i) of H° 
in energy order. Eigenstates of H, labeled by \a), are 
also ordered in energy, 

H\a)=E a \a), E a+1 > E a . (4) 

In order to obtain the GBWPT, let us divide the set 
of basis states into two parts, {\i), i = pi,pi + 1, • • - fa} 
and {\j}, j = 1, • • -pi — l,p-2 + 1, • • ■ N}. This also di- 
vides the Hilbert space into two sub-spaces, for which the 
corresponding projection operators are 

P2 

P=J2\i)(i\ and Q = 1- P. (5) 

i=Pl 

Subspaces related to the projection operators P and Q 
will be called in the following the P and Q subspaces , re- 
spectively. For an arbitrary eigenstate |a), which is split 
into two orthogonal parts \t) = P\a) and |/) = Q\a) by 
the projection operators P and Q, by making use of the 
stationary Schrodinger equation, it can be shown that 

I") = I*) + Ea l H o Q XV \ a )- ( 6 ) 
Introducing an operator T, 

T ^E-hp QXV > w 
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the iterative expansion of Eq. (^) is of the form 

\a) = \t) + T\t) + T 2 \t} + --- + T n - X \t) + T n \a). (8) 

Then, one can see that if the projection operators P and 
Q are such chosen that 



lim (a\(Tt) n T n \a) = 0, 



Eq. (H) gives 



|a) = \t) + T\t) + T 2 \t) +■■■ + T n \t) + 



(9) 



(10) 



Here the eigenvalue E a has been treated as a constant. 
Eq. (|lo| ) is a generalization of the so-called Brillouin- 
Wigner perturbation expansion (GBWPE) (for BWPE, 
see, for example, |p7[p8| ). 

Equation (|) is the condition for the GBWPE © to 
hold. Generally, if P and Q subspaces are such chosen 
that \E a — E°\ is large enough compared with bXV for 
any basis state \j) in the Q subspace, |T™ > will van- 
ish when n — > oo. Therefore, generally there are many 
choices for the P and Q operators ensuring Eq. ([[(j]) to 
hold. Among them, the projection operator(s) P with 
the minimum number of basis states and the correspond- 
ing operator(s) Q with the maximum number of basis 
states will be denoted by P a and Q a , respectively, and 
the corresponding \t) and |/) parts of the state \a) will 
be denoted by \t a ) = P a \a) and |/ Q ) = Q a \a), respec- 
tively. The \t a ) part of the state \a) will be called the 
non-perturbative (NPT) part of |a), and the \f a ) part 
called the perturbative (PT) part of \a). Correspond- 
ingly, the eigenfunction of the state |a), i.e., its compo- 
nents in the basis states, can also be divided into NPT 



is that the corresponding XU operator in the Q subspace 
does not have any eigenvector \v) with \Xu v \ > 1. 

As an application of the above results to the WBRM, 
let us study the component of an eigenstate |a) on a 
basis state \j) in the Q a subspace, C a j = (j\a), i.e., 
a component of the perturbative part |/ Q ). In what 
follows, for convenience, we use \j) to indicate a basis 
state in the Q a subspace, \i) to indicate a basis state 
in the P a subspace and \k) for the general case. Sup- 
pose j is larger than p2 for the P a subspace, specifically, 
Pi + mb < j < p2 + [m +1)6 with m being an inte- 
ger greater than or equal to zero. Noticing that a result 
of the band structure of the Hamiltonian matrix of the 
WBRM is that (j\(QV) n \t a ) is zero when n is less than 
(m+ 1), from Eqs. ( |l(i| ) and ( pi] ) we have 

C aj = (3\(T m+1 +T m+2 + ---)\t a ) 
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E n — H 



a(\U) m (l + \U + ---)Q\V\t a ). (14) 



Since QXV\t a ) is a vector in the Q a subspace, it can be 
expanded in the eigenstates of U, 



Q\V\t a )=Y,huW). 

V 

A result of this expansion and Eq. (|l3| ) is 



(15) 



(1 + XU + (At/) 2 + ■ • -)QXV\t a ) 
= ^2 hu{l + Xu v + {Xu u f + ■ ■ 



Xu v 



(16) 



and PT parts. Eq. (10) tells that \ f a ) can be expressed in 

terms of \t a ), E a , XV and H°. In what follows, generally Then, substituting Eq. Jl6| ) into Eq. (|l4"|), one has 



the GBWPE (10) will be used to discuss the perturbative 
parts of eigenstates only. 

To achieve a more explicit condition for Eq. (^) to hold, 
let us write T n in the following form, 



Cai - E r 



^£[T^o>>K^r. (it) 



J V 



rpfl 



where 



E n - H° 



U = QV 



(AC/)™- 1 • QXV, 



1 



E a — H 



(ii) 



(12) 



is an operator in the Q subspace. Eq. (JllJ) shows that it 
is the properties of XU that determines if T n \a) vanishes 
when n approaches to oo. Indeed, writing eigenstates of 
the operator U in the Q subspace as \v), 



U\v) = u v \u), 



(13) 



it is easy to see that if all the values of | Xu v \ are less than 
1, then T n \a) vanishes when n goes to infinity. Therefore, 
generally the condition for Eq. (|l0|) (also Eq. (0)) to hold 



Since \Xu v \ < 1 for all u, the behavior of the long tails 
of the EFs of the WBRM is more or less like exponential 
decay as discussed in Ref. Q. In fact, from the viewpoint 
of the GBWPT, the proof and arguments given in Ref. 
§ for behaviors of long tails of the LDOS of the WBRM 
are still valid when perturbation is strong. 

What is of special interest here is the behavior of C a j 
for j in the regions of {P21P2 + b] and [pi —b,p{), i.e., for 
j in the regions of size b just beside the non-perturbative 
part of the eigenstate. For these j, m = in Eq. jl^ ) 
and it is not necessary for \C a j\ to decay exponentially. 
That is to say, the decaying speed of \C a j\ for these j 
should be slower than that for j larger than P2 + b or 
less than p\ — b. The two regions of j 6 {p2,P2 + b] and 
j G [pi — b,pi) will be called the slope regions of the 
eigenstate \a) ( the reason for such a name will be clear 
from numerical results in the next section). 
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According to Eq. (|io|), the explicit expression for C a j 
in terms of elements of the Hamiltonian matrix can be 
written as 



P2 



XVjk XVki 



i=pi J feeQ J K 

XVjk X y kl XVu 
+ ^ Q E a — E" E a — El E a -E? ^ (18j 

where C a i = (i\t a ). This expression can be written in 
a simpler form by making use of the concept of path, 
in analogy to that in the Feynman's path integral theory 
p9|] , which can be done as follows. For (q+1) basis states 
|fci), • • • |fc g -i), \i)} with only \i) in the P a subspace, 



we term the sequence j — > k\ 



Cg-l 



i a path 



of q paces from j to i, if Vkk' corresponding to each pace 
is non-zero. Clearly paths from j to i are determined by 
the structure of the Hamiltonian matrix in H° represen- 
tation. Attributing a factor Vkk' l(E a — E^) to each pace 
k — > fc', the contribution of a path s to C a j, denoted by 
fasij — > *)i i s defined as the product of the factors of all 
its paces. Then, C a j in Eq. (fL8|) can be rewritten as 



/'2 



(19) 



An advantage of expressing C a j in terms of contribu- 
tions of paths is that it makes it easier to understand the 
important role played by the size of the non-perturbative 
part of the EF in determining its shape. Such a size for 
the state \a) will be denoted by N p , N p = p% — p± + 1. 
(N p is, in fact, a function of a, but, for brevity, the sub- 
script a will be omitted.) Two values of A are of special 
interest when A increases from zero. One is the smallest 
A for N p = 2, indicating the beginning of the invalidity of 
the ordinary Brillouin-Wigner perturbation theory. The 
other A is for the case of N p = b. The value of this A is 
of interest since topologically the structure of paths for 
the case of N p > b is different from that for N p < b. 
For example, when N p > 6, paths starting from j < p\ 
can never reach a j 1 which is greater than p2- That is 
to say, no path can cross the non-perturbative region. 
On the other hand, when N p < b, paths can cross the 
non-perturbative region. One can expect that such a dif- 
ference should have some effects on the shapes of both 
EFs and LDOS. 



III. NUMERICAL STUDY FOR THE 
PERTURB ATIVE AND NON-PERTURBATIVE 
PARTS OF EIGENFUNCTIONS 

In this section we study numerically the division of 
EFs into perturbative (PT) and non-perturbative (NPT) 
parts. Both the individual shape and the average shape 
of EFs will be studied. 



The boundary of the NPT part of a state \a), i.e., the 
values of pi and P2, are calculated in the following way. 
For a Hamiltonian matrix as in Eq. (^), we first diago- 
nalizc it numerically and get all its EFs. Then, for an 
eigenstate \a) with energy E a , we find out all the pairs 
of pi and P2 for which the value of (a\(T^) n T n \a) could 
become smaller than a small quantity (say 10 -6 ) when 
n goes to infinity. Finally, for the pair of (jpi,p2) thus ob- 
tained with the smallest value of N p = p2 — p\ + 1 and the 
pairs of (pi + 1 , p2 ) and (px , P2 — 1) , we calculate eigenval- 
ues of the corresponding U operators in Eq. ([l2]) in order 
to check out if all the \Xu v \ for the pair of (px,P2) are less 
than 1 while some of the (Am^I for each of the other two 
pairs (pi + 1,7*2) and (px,P2 — 1) are larger than 1. 

The shape of an eigenstate \a) in the unperturbed 
states can be defined as 



W a (E°)= y\C ak \ 2 5(E° - E° k ). 



k 



(20) 



In order to obtain the average shape of eigenstates, dif- 
ferent individual distributions W a (E°) should be first 
shifted into a common region. For this, we express 
W a (E°) with respect to E a before averaging. For a from 
oei to a 2 , with W a (E°) expressed as W a (E a — E a ), the 
average shape of the eigenstates, denoted by W(E°), can 
be calculated. Numerically, W(E®) are thus calculated: 
For a state |or), some value of E® and all k satisfying 
\E® — (El — E a )\ < 5E w /2, where SE W is some quantity 
chosen for dividing E® into small regions, we calculate 
the sum of J^k |Cafc| 2 i then take the average of the sum 
over a from ot\ to a.2- The result thus obtained is the av- 
erage value of W in the region (E° - SE w /2, E° + 5E w /2) 
and is also denoted by W(E®). Similarly the values 
of W(E° ± mSE w ) (m = 1,2,---) can also be calcu- 
lated. The average shape of eigenstates can also be di- 
vided into a NPT part and a PT part by the averaged 
boundary of the NPT parts of the states \a), denoted 
by Px = (pi - E a ) and p^ = (p 2 - E a ). The size 
of the NPT part of the average shape of eigenstates is 
(N p ) = (j>2 — Pi + 1). The values of A for (N p ) beginning 
to be larger than 1 and for (N p ) = b will be denoted by 
A / and Af, , respectively. 

Now let us present the numerical results. The first one 
is for the case of N = 300, b = 10 and A = 0.6. This is 
a case for which N p can be equal to both 1 and 2. The 
values of |C Q fe| 2 = |(/c|o;)| 2 for four EFs in the middle 
energy region (a = 148 — 151) are given in Fig. 1. The 
two vertical dashed-dot straight lines for each case indi- 
cate the positions of pi and P2 of the NPT part of the 
eigenstate, i.e., the boundary of the NPT part. Since 
Efc = k, pi and P2 are also the unperturbed eigenen- 
ergies of the corresponding basis states \k = pi) and 
\k = P2)- For a = 148 and 149, since p\ — P2 = a 
there is only one vertical dashed-dot line in each case. 
For a — 150 and 151, there are, in fact, two NPT parts 
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for each case. In Fig. 1 we give one of them for each case 
only, i.e., (pi,P2) = (149, 150) for a = 150 (the other one 
is (150,151)) and (pi,p 2 ) = (151,152) for a = 151 (the 
other one is (150,151)). Numerically we have found that 
only for small A is it possible for an eigenstate to have 
more than one NPT parts. The average shape of EFs for 
a from 130 to 170 is given in Fig. 2 with the correspond- 
ing values of pi and p^ indicated by vertical dashed-dot 
lines. SE W = 1 for calculating this W(E®) function. 

As mentioned in section II, another perturbation 
strength of interest is for the case of N p being able to 
equal the band width b of the Hamiltonian matrix. An 
example for this case is A = 1.4. For this value of 
A, there is only one NPT part for each state |a), e.g., 
N p = 9,10,7,9 for a = 148,149,150,151, respectively. 
Individual EFs for a = 148 — 151 are presented in Fig. 3 
with their boundaries of the NPT parts. Although for 
an individual EF the value of |C a fc| 2 could be still large 
outside the NPT region (Fig. 3), the main body of the 
average shape of EFs lies obviously in the averaged NPT 
region as shown in Fig. 4(a). For this averaged EF we 
still have 8E W equal to 1 and a from 130 to 170. An 
interesting feature of the distribution InW in the slope 
regions [p\ — b,p1) and (p%,p% + b] in Fig. 4(b) is that, 
as predicted in section II, the decaying speed in the two 
slope regions is obviously slower than that in the long 
tail regions. From Fig. 4(a) it is quite clear why we call 
the regions "slope" . 

In Ref. |p| the average shape of EFs for different val- 
ues of the Wigner parameter q, 



E a should be close to the centers of the corresponding 
NPT regions, that is, 



{pvf 



(21) 



and an ergodicity parameter "A e " (denoted by A in that 
paper), 



with aw 1.2, 0.92, 



(22) 



has been studied numerically. Results in that paper show 
that when q is small and A e not small, e.g., q = 1 and 
A e = 3.7, two averaging procedures, namely, averaging 
with respect to centers of energy shell, which are close to 
the perturbed eigenenergies, and averaging with respect 
to centroids of EFs, give similar results for the average 
shape of EFs; while for large q and small A e , e.g., q — 90 
and A e = 0.24, the two averaging procedures give differ- 
ent results and a so-called localization in energy shell has 
been found. 

What is of interest here is to study the difference be- 
tween the above two cases of different parameters q and 
A e from the viewpoints of the GBWPT. In Fig. 5, four 
examples of individual EFs are given for A = 4.0, for 
which q = 1.6 and A e = 1.82. When edge effects can be 
neglected, from the condition for Eq. (^) to hold given in 
the above section, it can be shown that the eigenenergies 



E n 



P2 - Pi 



(23) 



Here by the NPT region of an EF we mean the region of 
E^ between the boundaries (dashed-dot lines) of the NPT 
part of the EF as shown in Fig. 5. Then, since as shown 
in Fig. 5 centroids of the EFs are generally not far from 
the centers of the NPT regions, they are generally not far 
from the corresponding eigenenergies. Therefore, in this 
case the averaging with respect to the centroids of the 
EFs and the averaging with respect to the eigenenergies 
do not have much difference, and the average shape of 
EFs obtained by the two methods should be more or less 
similar. However, since there are still some differences 
between the centroids of the EFs and the eigenenergies, 
results of the two averaging methods can not be quite 
close to each other. In fact, as shown in Ref. results 
of the method of averaging with respect to centroids of 
EFs should be a little narrower than those of the other 
method. 

The average shape of EFs for a from 130 to 170 and 
A = 4.0 is given in Fig. 6 with SE W = 3.0. (In this pa- 
per numerically we study the averaging with respect to 
eigenenergies only.) For this averaged EF the slope re- 
gions are even clearer. The sizes of the two slope regions 
in Fig. 6 are larger than b, the predicted one for a sin- 
gle EF, since the values of N p for the NPT parts of the 
eigenstates \a) taken for averaging are not the same. In 
fact, the variance of the N p for the states is about 1.56. 

For studying the case of localization in energy shell as 
in Ref. Jp|], the dimension N = 300 is too small. There- 
fore, we increase the value of N to 900 and study the 
case of b = 10 and A = 30.0, for which the parameters 
q and A e are equal to 90 and 0.24, respectively, as those 
in Ref. ||. The values of \C ak \ 2 (dots) of four EFs in 
the middle energy region are given in Fig. 7 with the cor- 
responding boundaries (p± and P2) of the NPT parts of 
the EFs (indicated by vertical dashed-dot lines) . Clearly 
the main body of each of the EFs occupies only part of 
the NPT region. Then, resorting to Eq. (p3|), it is easy 
to see that for the average shape of EFs the result of 
averaging with respect to centroids of EFs should be ob- 
viously narrower than that of averaging with respect to 
eigenenergies. This phenomenon is called localization in 
energy shell in Ref. p3| . In our opinion, it is, in fact, 
localization of EFs in their NPT regions. 

Although, as shown in Fig. 7, when A = 30.0 many 
components of the EFs are quite small in the NPT parts, 
they can still be distinguished from those in the PT parts. 
In Fig. 8 we present the values of |C a fc| 2 in Fig. 7 in loga- 
rithm scale, from which the difference between the \C a k | 2 
in the NPT parts and those in the PT parts is quite 
clear. In fact, it can be seen from Fig. 8 that |C Q fc| 2 
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for small components in the NPT parts also decay ex- 
ponentially on average, but the decaying speed is slower 
than that for components in the PT parts. Similar re- 
sults have also been found for other eigenfunctions and 
the centroids of the EFs have been found scattered in 
the NPT regions. Therefore, averaging with respect to 
eigenenergies (around the centers of the NPT regions) is 
reasonable when studying the average shape of EFs. The 
average shape of EFs in the middle energy region for a 
from 420 to 480 is shown in Fig. 9 (8E W = 10). The 
two slope regions are also obvious. The widths of the 
two slope regions, which can be seen from the figure is 
about 28E W = 2b, are also larger than that predicted for 
a single EF due to the fact that the variance of N p for 
these states \a) is about 5.14. 

In order to have a clear picture for the variation of the 
average size (N p ) of NPT parts of EFs with the perturba- 
tion strength A, we plot it in Fig. 10(a) (for N = 300 and 
a from 130 to 170). The two solid straight lines in the 
figure show that (N p ) has a good linear dependence on A 
in the two regions of A £ (2.0, 4.5) and A e (4.5, 7.0). The 
variation of (N p ) for A less than 2 is given in the inset of 
the figure in more detail. The solid curve in the inset is 
a fitting curve of a quadratic form: (7.5(A — 0.4) 2 + 1.0). 
Since in this case the value of A/ is a little larger than 0.4, 
i.e., for A < 0.4 the size N p of the NPT parts of the EFs 
is equal to one, the region of A < 0.4 has been neglected 
when fitting the (N p ) by the quadratic curve. From the 
inset it can be seen that the fitting curve fits the values 
of (N p ) quite well in the region of A £ (0.4, 1.5). Inter- 
estingly, the value of A& for (N p ) — bis between 1.45 and 
1.5, that is to say, in the region of A £ (A/, A&), the values 
of (N p ) have a quadratic dependence on A. 

The variance of N p , denoted by AN p , is given in 
Fig. 10(b) (circles), which shows that AN P is small com- 
pared with N p . From the figure it can be seen that for 
A < 5 the value of AN p increases with A on average. 
Then, since, as mentioned above, the value of AN p is 
about 5.14 when A = 30, it does not change much in the 
region of A £ (5, 30), that is to say, the relative value of 
the variance, namely AN p /N p , becomes even smaller as 
A becomes larger. 



IV. RELATIONSHIP BETWEEN THE SHAPE OF 
EIGENFUNCTIONS AND OF LDOS 

The existence of a relationship between the average 
shape of EFs and that of LDOS in some particular cases 
has already been noticed numerically by several authors 
(see, for example, ||Jl|||). Here for the WBRM we 
show that some analytical arguments can be given for it 
by making use of the GBWPT. 

The so-called local spectral density of states (LDOS) 
for an unperturbed state |fc) is defined as 



p k L (E)=Y,\ C ^\ 2 KE-E a ) (24) 

a 

where C a k = (k\a). Making use of the division of EFs 
into NPT and PT parts, the non-perturbative (NPT) and 
perturbative (PT) part of the LDOS can be defined in the 
following way. Let us consider all the perturbed states 
| a) for which C a k is in the NPT parts of their eigenfunc- 
tions. Suppose a Pl is the smallest one among these a 
and the largest one among them, i.e., \a Pl ) has 

the smallest eigenenergy and \a P2 ) has the largest. A 
property of the WBRM is that for any state \a) with 
ot Pl < ol < Oi P2 , generally to say, C a k is in the NPT part 
of its eigenfunction. Then, the non-perturbative (NPT) 
part of the LDOS p\(E) can be defined as the sum in 
Eq. (HJ) over a £ [a Pl ,a P2 ], and the perturbative (PT) 
part can be defined as the sum over the rest a. The size 
of the NPT part of the LDOS is defined as (a P2 — a pi +1) 
and will be denoted by N a . 

The average shape of LDOS, denoted by pl(E s ), can 
be obtained in a way similar to that for the EFs discussed 
in section III, except that for the LDOS the distributions 
p\(E) are expressed as functions of (E — Efy before av- 
eraging. Since E® is also the centroid of p\{E), such an 
averaging method is the most natural one. The average 
shape of LDOS can also be divided into a NPT part and 
a PT part by the averaged boundary of the NPT parts 
of the individual LDOS. 

Before discussing the relation between the average 
shape of EFs and that of LDOS, let us first draw some 
conclusions for properties of EFs from the numerical re- 
sults discussed in section III. Firstly, based on the nu- 
merical results given in section III for the average shape 
of EFs, we can make such an approximate treatment for 
the average values of |C a i| 2 in the NPT parts of EFs: 
for a given A not extremely large, they can be treated 
as a constant denoted by t\, i.e., |C„j| 2 w i 2 , when edge 
effects can be neglected. Here we would like to stress 
that Figs. 2,4,6,8 show that the values of |C Q i| 2 on the 
edges of the NPT parts become smaller than those in the 
middle regions when A increases and the above approx- 
imate treatment may fail in case of extremely large A. 
Secondly, another result of the last section is that the 
variance AN p is small compared with N p . This means 
that when discussing approximate relations it is reason- 
able to treat the size N p of NPT parts of eigenstates as 
a constant for a fixed A. 

As results of the above two approximations and the 
approximate relation (|2^), it can be shown that, under 
the corresponding conditions, (a) the average value of 
\C ak \ 2 in the NPT parts of LDOS is also about t\, (b) 
the size of the NPT part of a LDOS is close to that of 
an EF, i.e., N ap w N p , (c) usually E® is in the mid- 
dle region of the NPT part of the LDOS p k L (E), i.e., 
El — E apl sa E ap2 — E®. Then, one can see that the 
following approximate relation holds for the NPT part 



G 



of an averaged EF W and the NPT part of an averaged 
LDOS PL , 

W(E) _ pl{EI 

P°(E) ~ p{E) [2b) 

where p°(E) and p(E) are the averaged density of states 
of the unperturbed and perturbed spectra in the corre- 
sponding regions, respectively,. Here, for brevity, instead 
of E® and E s , we use E to denote the variables of the 
functions W and p^. 

Next let us discuss the relation between the PT part of 
the average shape of EFs and the PT part of the average 
shape of LDOS. According to Eq. (|l9|), for a basis 
state \j) in the Q a subspace can be written as 

P2 P2 

^aj = fas 1 (j — > il)fas 2 {j ~ * ai\C <xi 2 - 

il=Pl »2=Pl S!,S 2 

(26) 

Following a way similar to that of calculating the average 
shape of EFs W(E S ) in the last section, we take the av- 
erage of this quantity over different eigenfunctions. Since 
C a i in the NPT parts of EFs have random signs, the main 
contribution to in Eq. ( p6|) comes from the i\ = i% 
terms and 

P2 

C lj * E E /a.i(j "►<)**■ (27) 

i=pi si,s 2 

When the number of eigenstates taken for averaging is 
large enough, since f as is a product of factors Vkk' / (E a — 
E^), the main contribution in ( |27j ) to C^j comes from the 
terms with path si equal to path S2, that is, 

i=pi s 

Let us first discuss the relation (^), which is for cases 
without enough eigenstates taken for averaging. In these 
cases, due to the signs of the denominators (E a — ET) 
of the factors of paces, there is a systematic difference 
between the contribution of paths starting from j with 
E® < E a and the contribution of paths starting from j' 
with E®, > E a , especially when N p > b. A result of this 

difference is that |C Q:) | 2 for the left PT part of an aver- 
aged EF may be different from the corresponding one for 
its right PT part. But, for PT parts of averaged EFs on 
the same side, such a difference does not exist. In fact, 
when edge effects can be neglected, the structure of paths 
is similar for different EFs when the difference in the size 
of NPT parts of the EFs can be neglected, and we have 

f asi (j -> i)f aS2 {j -> i) ~ fa's^W i')f*'s 2 {j' i') 

(29) 



where s[ s' 2 are paths similar to s\ and s 2 , respectively, 
but from j' to i' with f — i' = j — i and j' — p\ — j — p\. 
From the relations ( p7| ) and (p9|), one can see that, when 
edge effects can be neglected, the left (right) PT parts of 
averaged EFs in different energy regions should be simi- 
lar. By definition, the left (right) PT part of the average 
shape of EFs should be related to the right (left) PT 
part of the average shape of LDOS, since, for example, 
a \C aj \ 2 in the left PT part of the EF W a (E°) is in the 
right PT part of the LDOS p j L {E). Then, the relation- 
ship between the PT part of the average shape of EFs 
and that of the average shape of LDOS should be 

W(E) ^ p^ L (E) 

P°(E)~ P(E) {M) 

where P-l(E), the inverted LDOS, is defined as 
P-l(E) = pl(—E). This means that the left (right) 
tail of the average shape of EFs is more similar to the 
right (left) tail of the average shape of LDOS than its 
right (left) tail does. Since |C Q j| 2 can be approximately 
treated as a constant for the NPT parts of both EFs and 
LDOS, the approximate relation (^o|) also holds for the 
NPT parts of the average shape of EFs and of LDOS. 

When there are enough eigenstates taken for averag- 
ing, we have the relation (|2^). In this case, instead of 
(p9|), we have 



flU^i)*& s >U'^i'), (31) 

and the left PT part of an averaged EF should be similar 
to its right PT part. Then, the same relation as in ( |25| ) 
can be obtained for the PT part of the average shape of 
EFs and the PT part of the average shape of LDOS. In 
this case, pl(E s ) « p- L (E s ). 

The above results have been checked by numerical cal- 
culations. In Fig. 11(a) and (b), an average shape of 
EFs W(E) (circles) is compared with a corresponding 
inverted average shape of LDOS P-l{E) (solid curve) 
for A = 4.0 and b = 10. In order to avoid edge effects, 
N is chosen to be 900. In this case, in the middle en- 
ergy regions the difference between p(E) and p°(E) can 
be neglected. The averaging has been done for 60 per- 
turbed and 60 unperturbed states in the middle energy 
regions, respectively. For this relatively small number of 
states taken for averaging, one can expect that the re- 
lation (^5|) is not so good as the (J3TJ) in the tail region. 
Indeed, as shown in Fig. 11(c) the LDOS Pl(E) (solid 
curve) is not so close to the W(E) (circles) as the in- 
verted one P-l(E) in Fig. 11(b). In order to show the 
influence of edge effects on the relation (|30|), in Fig. 11(d) 
we give a result obtained when N is equal to 300. 

Finally, we would like to stress that although the above 
arguments leading to the approximate relations ( |25| ) and 
(p0|) are for not very large A, numerical results in Ref. 
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show that central parts of averaged EFs (for averag- 
ing with respect to eigenenergies) are still close to those of 
averaged LDOS when A is quite large. In fact, for a given 



A, the approximation \C a 



t{ for NPT parts of EFs, 



which has been employed in the above arguments, is not 
necessary for deducing the relations ( p5j) and ( |30|). For 



example, under the approximation that \C mi | 



I Cm 2 1 2 



with i\ —pi = P2 — *2 for NPT parts of EFs, similar argu- 
ments as given above also lead to the two relations ( |25| ) 
and (|30|). In fact, this approximation is in even better 
agreement with the numerical results given in the last 
section than the former approximation |C Q i| 2 ~ t\ and 
may be still valid for very large A. Therefore, the two 
relations ( p5|) and (|3^) may be correct for very large A, 
too. 



V. VARIATION OF THE SHAPE OF LDOS WITH 
PERTURBATION STRENGTH 

Analytical and numerical study for the shape of the 
LDOS of the WBRM has been done in, e.g., 
It is known that when perturbation is weak (but not very 
weak) the average shape of LDOS is of the Breit-Wigncr 
form (Lorcntzian distribution), and when perturbation 
becomes strong the shape will approach to a form pre- 
dicted by the semi-circle law. 

Our interest here is in studying the process of the tran- 
sition from the Breit-Wigner form to the semi-circle form, 
in order to see if the GBWPT can throw light on how the 
transition occurs. In fact, due to two results of the previ- 
ous two sections that (a) the central part of the average 
shape of EFs is composed of its NPT part and the slope 
region of the PT part and (b) the average shape of EFs is 
similar to that of LDOS when the density of states of the 
perturbed spectrum is similar to that of the unperturbed 
spectrum, one can expect that properties of the average 
shape of LDOS, particularly of its half-width, should be 
related to properties of (N p ), the average size of NPT 
parts of EFs. 

An interesting property of the semicircle law 



Psc{E) 



Rl 



E\ 



\E\ < R, 



o- 



(32) 



where Rq = A V 8b, is that it obeys a scaling law. Specif- 
ically, after a rescaling 

R = XR'o, E = XE', Psc (E) = p' sc (E')/X, (33) 

it becomes 

2 



P'sc(E') 



ttR 



. , y R' Q E' 



(34) 



having the same form as the p sc in Eq. (|32[). In fact, 
Eq. (34) is just Eq. for the case of A = 1. This prop- 
erty of the semicircle law supplies a convenient method of 



studying the approaching of the average shape of LDOS 
Pl{E s ) (for brevity, as in the last section, the subscript s 
for the variable E s will be omitted in what follows) to the 
semicircle form p sc (E) with increasing A. That is, first 
we change the LDOS Pl{E) for a perturbation strength 
A to a rescaled one pY(E) defined by 



pl s (E) = X PL (XE), 



(35) 



then, we compare it with the semicircle form p sc (E) for 
A = 1. 

In Fig. 12 we present such comparisons for A from 1.2 
to 1.9. The LDOS Pl(E) in this and the following figures 
are for N = 500 and b = 10. In order to take average, 
30 Hamiltonian matrices of different realizations of the 
random numbers have been diagonalized for each A and 
50 individual LDOS p k L {E) in the middle energy region 
have been taken for each Hamiltonian matrix. Figure 
12 shows that the sign for the p r ^{E) to approach the 
semicircle form appears when A is between 1.5 and 1.6. 
Interestingly, in this case the value of A& for (N p ) = b is 
also between 1.5 and 1.6. As indicated in section II, such 
a value of A is of interest, since for a perturbed eigen- 
state with N p > b paths starting from the left PT part 
of the state can not reach the right PT part, and vice 
versa. That is to say, when A becomes larger then Xb, 
there will be a change in the topological structure of the 
paths. Therefore, such a coincidence should not be an 
accident. 

When A becomes larger, as is known, the form of 
pY(E) will become closer to the semicircle form (Fig. 13). 
Since the semicircle law has a scaling behavior, an in- 
teresting result of the comparisons between the rescaled 
LDOS p r ^{E) and the semicircle form given in Figs. 12 
and 13 is that the LDOS Pl(E) obeys an approximate 
scaling law when A is larger than At>. As a consequence, 
the dependence of the half-width of the LDOS on A 
should become linear when A exceeds A;,. 

The Brcit-Wigner form 



Pbw(E) 



r/27r 



e 2 + r 2 /4 



(36) 



does not have the scaling property as the semicircle law. 
For the purpose of studying the relationship between the 
Breit-Wigner form and the average shape of LDOS, we 
use the former as a fitting curve for the later with the 
width r as the fitting parameter. The fitting is done by 
requiring the minimum of the area difference 



AS = J \p L (E) - PBW {E)\dE 



(37) 



between the Breit-Wigner form and the histogram of 
the averaged LDOS. In numerical calculations, we first 
change the Breit-Wigner form to a histogram correspond- 
ing to the histogram of the LDOS, then calculate the AS. 
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Four examples thus obtained for A = 0.4,0.7,1.0 and 
1.5 are given in Figs. 14 and 15. From numerical results 
we have found that the average shape of LDOS begins 
to be fitted well by the Breit-Wigner form just before A 
reaches 0.4. Interestingly, the value of A/ for (N p ) be- 
ginning to be larger than 1 is a little larger than 0.4. 
That is to say, the LDOS of the Breit-Wigner form be- 
gins to appear just before the ordinary Brillouin-Wigncr 
perturbation theory fails. With the increasing of A, the 
closeness between the LDOS and the Breit-Wigner form 
maintains for a small region of the A. Then, when A in- 
creases further, deviation will become more obvious. In 
fact, when A = A& 1.5, the LDOS Pl(E) is absolutely 
different from the Breit-Wigner form, while it becomes 
closer to the semicircle form. 

Since the smallest area difference AS in Eq. (^) gives 
a measure for the deviation of the average shape of 
LDOS from the Breit-Wigner form, we plot it in Fig. 16 
(squares). It can be seen that the deviation reaches its 
saturated value at about A = 2.0. Similarly, in order to 
show the deviation of the average shape of LDOS from 
the semicircle form, one can introduce another area dif- 
ference 

A5 = J \ P r L s (E)- Psc (E)\dE (38) 

where p sc (E) is for A = 1. Variation of this AS" with 
A is also given in Fig. 16 (circles). It shows that before 
A reaches 1.5 the deviation is large and drops quickly. 
For A from 1.5 to about 4.5, AS drops slower. When A 
is larger than 4.5, the values of AS* are quite small and 
change quite slowly, which means that the average shape 
of LDOS has become quite close to the semicircle form. 

Variation of the half-width of the average shape of 
LDOS is of both experimental and theoretical interest. 
Such a variation is given in Fig. 17 by triangles. As ex- 
pected, when A exceeds Xb ~ 1.5, the half-width can be 
fitted well by a straight line. The variation of (N p ) with 
A is also given in Fig. 17 (circles). Corresponding to the 
three regions of A for the variation of the AS measur- 
ing the deviation of the LDOS from the semicircle form 
represented by circles in Fig. 16, the variation of (N p ) 
can also be divided into three regions: (1) A < 1.5, (2) 
1.5 < A < 4.5, in which it can be fitted well by a straight 
line, and (3) A > 4.5, in which it can be fitted well by 
another straight line. 

The fitting for the half-width of the LDOS and for the 
(N p ) for small A is given in the upper-left inset in Fig. 17. 
Since when A < 0.4 the value of (N p ) equals to 1, the fit- 
ting for (N p ) (circles) has been done for A > 0.4 by a 
quadratic curve (a(A - 0.4) 2 - 1.0) with a w 7.0. The 
quadratic feature of the dependence of (N p ) on A is quite 
clear in the region of 0.4 < A < 1.5. For the half-width 
(triangles), the same form of fitting curve with a « 15.0 
is also given in the inset. As expected, the quadratic 
feature is also clear for the half-width in the region of 



0.4 < A < 1.5. Therefore, the correspondence between 
behaviors of the average size of NPT parts of EFs and 
those of the half width of the average shape of LDOS is 
quite clear. The lower-right inset in Fig. 17 gives a com- 
parison between the (N p ) for TV = 500 (solid curve) here 
and the (N p ) for N = 300 (dashed line) in Fig. 10. The 
small deviation between them comes from both edge ef- 
fects and the fact that they are for Hamiltonian matrices 
with different values of off-diagonal elements. 

In conclusion, as expected, numerical results given in 
this section show that properties of the average shape of 
LDOS are indeed related to properties of the average size 
of NPT parts of perturbed eigenstates and knowledge of 
the latter, especially the value of A& for (N p ) = b, can 
indeed give deeper understanding for the former. 

VI. CONCLUSIONS 

The Wigner Band Random Matrix (WBRM) model is 
studied numerically in this paper by making use of a gen- 
eralization of the Brillouin- Wigner perturbation theory 
(GBWPT). According to the GBWPT, an energy eigen- 
function (EF) of a perturbed system can be divided into 
a non-perturbative (NPT) and a perturbative (PT) part 
with the PT part expressed as a perturbation expansion. 
Further more, the PT part can be divided into a slope 
region and a tail region. 

Numerically we have found that for the average shape 
of EFs its central part is composed of its NPT part and 
the slope region of its PT part. That is, the GBWPT 
can give an analytical definition for the division of the 
average shape of EFs into central parts and tails. For 
the shape of individual EFs, numerical results show that 
when the perturbation is not strong their NPT part and 
the slope region of their PT part are usually composed of 
large components. But when the perturbation becomes 
stronger, the region of the NPT part of an EF occupied 
by large components will become relatively smaller, i.e., 
the ratio of the region to the whole NPT part will be- 
come smaller. As for the small components in the NPT 
part of an EF in this case, the difference between them 
and the components in the PT part is also obvious. Here 
we would like to point out that the possibility of dividing 
EFs into NPT and PT parts should be useful in reduc- 
ing calculation time for eigenfunctions in a given energy 
region. 

It is already known from some previous numerical stud- 
ies that there is a relationship between the average shape 
of EFs and that of LDOS, but the reason is not clear. 
Resorting to the GBWPT and some conjectures made 
from numerical results, it is possible to show that such 
a relationship indeed exist for the WBRM. Particularly, 
when the number of states taken for averaging is not large 
enough, the relationship is between the average shape of 
EFs and the average shape of inverted LDOS. 
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A result of the above properties of the average shape 
of EFs and of LDOS is that some properties of the aver- 
age shape of LDOS should be related to properties of the 
NPT part of the average shape of EFs. This has been 
studied in detail by numerical calculations. Firstly, it is 
found that the LDOS of the Breit-Wigner form appears 
just before the perturbation strength A reaches a value 
A/, for which the NPT parts of EFs begin to have more 
that one components. Secondly, when A reaches At,, for 
which the average size (N p ) of NPT parts of EFs is equal 
to the band width b of the Hamiltonian matrix, the av- 
erage shape of LDOS begins to be close to the semicircle 
form predicted by the semicircle law, and for A larger 
than Ab the average shape of LDOS obeys an approxi- 
mate scaling law. Thirdly, variation of the half-width of 
the average shape of LDOS with perturbation strength 
is closely related to the variation of the average size (N p ) 
of NPT parts of EFs. Particularly, when A < Ab, i.e., 
(N p ) < 6, both of them are of quadratic form, and when 
A is larger than Ab, both of them becomes linear. 
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FIG. 1. Circles show the values of \C a k | 2 = | (k\a) | 2 for four 
energy eigenstates \a) when A = 0.6 and N = 300. The verti- 
cal dashed-dot lines indicate the boundaries p\ and P2 of the 
corresponding four non-perturbative parts of the eigenstates 
(for a = 148 and 149, pi = p2 — ct). 
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FIG. 2. (a) The average shape of eigenstates in the middle 
energy region, W(E 3 ) (circles), for A = 0.6 and N = 300. 
The vertical dashed-dot lines indicate the averaged boundary 
Pi and P2 of the corresponding non-perturbative parts of the 
eigenstates. (b) Same as in (a) in logarithm scale. 



FIG. 3. Same as in Fig. 1 for A = 1.4. 



FIG. 4. Same as in Fig. 2 for A = 1.4. 



FIG. 5. Same as in Fig. 1 for A = 4.0. 



FIG. 6. Same as in Fig. 2 for A = 4.0. 



FIG. 7. Same as in Fig. 1 for A = 30.0 and N = 900. 
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FIG. 8. Same as in Fig. 7 in logarithm scale. 



FIG. 9. Same as in Fig. 2 for A = 30.0 and N = 900. 

FIG. 10. (a) Variation of the average size (N p ) (circles) of 
non-perturbative parts of eigenstates in the middle energy re- 
gion with the perturbation strength A for N = 300. The two 
solid straight lines are fitting lines. The inset shows the fitting 
curve of the form 7.5(A - 0.4) 2 + 1.0 for A from 0.4 to 1.5. (b) 
Variation of the variance of N p , AN P (circles), with A. 

FIG. 11. (a) A comparison between an average shape of 
eigenfunctions W(E) (circles) and a related average shape of 
inverted LDOS P-l(E) (solid curve) for A = 4.0 and N = 900. 
(b) Same as in (a) in logarithm scale, (c) A comparison be- 
tween the W(E) (circles) in (a) and the LDOS Pl(E) (solid 
curve) related to the P-l(E) in (a) in logarithm scale, (d) 
Same as in (b) for N = 300. 

FIG. 12. Comparisons between the rescaled LDOS p r L s 
(histograms) obtained when N = 500 and the prediction of 
semicircle law (solid curve) for A = 1. 



FIG. 13. Same as in Fig. 12 for larger values of A. 

FIG. 14. The average shape of LDOS Pl(E) (histograms) 
obtained when N = 500 and the fitting curves (dashed curves) 
of the Breit-Wigner form. 



FIG. 15. Same as in Fig. 14 for larger values of A. 

FIG. 16. Variation of the area difference AS (circles) be- 
tween the rescaled LDOS and the semicircle form for A = 1, 
and variation of the area difference AS (squares) between the 
average shape of LDOS and the fitting Breit-Wigner curves 
with A for N — 500. 

FIG. 17. The triangles show the values of the half-width of 
the average shape of LDOS, and the circles show the values of 
the average size (N p ) of non-perturbative parts of eigenfunc- 
tions when N = 500. The three solid straight lines are fitting 
lines. The upper-left inset shows the fitting curve of the form 
(7(A - 0.4) 2 + 1.0) for (N p ) and the fitting curve of the form 
(15(A - 0.4) 2 + 1.0) for the half-width of the LDOS for A from 
0.4 to 1.5. The lower-right inset gives a comparison between 
the (N p ) in Fig. 10 obtained when N = 300 (dashed curve) 
and the {N p ) here obtained when N = 500 (solid curve). 
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